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ABSTRACT 

A (sub-)millimeter line and continuum study of the class I protostar Elias 29 in the p Ophiuchi molec- 
ular cloud is presented, whose goals are to understand the nature of this source, and to locate the ices 
that are abundantly present along this line of sight. Within 15-60" beams, several different compo- 
nents contribute to the line emission. Two different foreground clouds are detected, an envelope/disk 
system and a dense ridge of HCO+-rich material. The latter two components are spatially separated 
in millimeter interferometer maps. We analyze the envelope/disk system by using inside-out collapse 
and flared disk models. The disk is in a relatively face-on orientation (< 60°), which explains many 
of the remarkable observational features of Elias 29, such as its flat SED, its brightness in the near 
infrared, the extended components found in speckle interferometry observations, and its high velocity 
molecular outflow. It cannot account for the ices seen along the line of sight, however. A small fraction 
of the ices is present in a (remnant) envelope of mass 0.12-0.33 Mq, but most of the ices ('^70%) are 
present in cool (T <40 K) quiescent foreground clouds. This explains the observed absence of thermally 
processed ices (crystallized H2O) toward Elias 29. Nevertheless, the temperatures could be sufhciently 
high to account for the low abundance of apolar (CO, N2, O2) ices. This work shows that it is crucial 
to obtain spectrally and spatially resolved information from single-dish and interferometric molecular 
gas observations in order to determine the nature of protostars and to interpret infrared ISO satellite 
observations of ices and silicates along a pencil beam. 

Subject headings: dust, extinction — Infrared: ISM — ISM: molecules — stars: formation — stars: 
individual (Elias 29) — submillimeter 



1. INTRODUCTION 

A rich chemical and physical interplay exists between 
gas and grains in which molecules are formed on grains, 
creating ice mantles that are preserved in environments 
ranging from quiescent dense molecular clouds to en- 
velopes and disks around protostars. Various processes, 
among which are bombardment by cosmic rays, ultraviolet 
irradiation, heating, and shocks, can physically or chem- 
ically alter the icy mantles, or return molecules into the 
gas phase (see Tielens & Charnley 1997; van Dishoeck & 
Blake 1998, and references therein). 

A study of the chemical evolution of dense clouds to 
planet forming disks would ideally involve observations of 
molecular gas and ices in a range of environments, from 
quiescent clouds to disk dominated protostars. The most 
pristine, initial conditions are presumably well sampled 
by field stars behind clouds tracing quiescent molecular 
cloud material (e.g. Whittet et al. 1998). Lines of sight 
to protostars are more difficult to characterize, however, 
since they may trace quiescent foreground material, in ad- 



dition to the gas and ices in their envelopes and disks (e.g. 
Boogert et al. 2000b). It is thus crucial to character- 
ize the line of sight conditions in order to locate the ices 
and derive physical conditions and eventually the evolu- 
tion of the molecular gas and solid state in the interstellar 
medium. 

In this paper, we study the line of sight of the class I 
protostar Elias 29 in the p Oph molecular cloud, using 
(sub)millimeter single dish and interferometer gas phase 
observations. This object is one of the most luminous pro- 
tostars (~ 36 Lq\ Chen et al. 1995) in the nearby p Oph 
complex (d ~ 160 pc; Whittet 1974), yet little is known 
about its nature and line of sight conditions. Abundant 
ice has been detected in its direction (Zinnecker, Webster, 
& Geballe 1985; Tanaka et al. 1990). A detailed analysis 
of ice band profiles indicates that the ices are not strongly 
thermally processed (i.e. the ices are not crystallized or 
segregated) , despite the presence of abundant warm molec- 
ular gas toward the object (Boogert et al. 2000b). This 
contrasts with high mass, luminous (> 10^ L0) proto- 
stars, where significant thermal processing of the ices ac- 



^ California Institute of Technology, Downs Laboratory of Physics 320-47, Pasadena, CA 91125, USA; boogertOsubmm. caltech.edu 
2 RAL, University of California at Berkeley, Astronomy Department, 601 Campbell Hall # 3411, Berkeley, CA 94720, USA 
current address: Steward Observatory, University of Arizona, 933 N. Cherry Ave., Tucson, AZ 85721, USA 

* Observatoire de Bordeaux, B.P. 89, 33270 Floirac, France 

^ Kapteyn Astronomical Institute, P.O. Box 800, 9700 AV Groningen, the Netherlands 

^ SRON, P.O. Box 800, 9700 AV Groningen, the Netherlands 

^ Leiden Observatory, P. O. Box 9513, 2300 RA Leiden, the Netherlands 

* California Institute of Technology, Division of Geological and Planetary Sciences 150-21, Pasadena, CA 91125, USA 
9 California Institute of Technology, SIRTF Science Center, IPAC, Pasadena, CA 91125, USA 



2 



Boogert et al. 



Table 1 
Observational Summary 



Telescope Beam (") \ oi f Species" Date or Reference 



NRAO-12m 


43-86 


70-145 GHz 


CO, CS, HCO+ 


05/1995 








H2CO, CH3OH 


JCMT 


14-22 


200-400 GHz 


CO, CS, HCO+ 


1995-1997 








H2CO, CH3OH 




JCMT 


7 


692 GHz 


CO 6-5 


Ccccarclli et al., in prep. 


CSO 


21-35 


200-400 GHz 


H2CO, CH3OH 


1999-2001 








CO, HCO+ maps 




OVRO 


4x8, 3x6 


87, 110 GHz 


CO, HCO+, SiO 


09/1999-07/2000 


OVRO 


4x8, 3x6 


2.7, 3.3 mm 


Continuum 


09/1999-07/2000 


IRAM-30m 


15 


1.3 mm 


Continuum 


Motte et al. 1998 


ISO SWS 


14-33 


2.3-45 nm 


Continuum 


Boogert et al. 2000b 








Ices. Silicates 




ISO LW S 


80 


15 200 //m 


( Vjiil iiiuuiii 


LSooscrt rl al. 20001) 



^ main species and/or isotopes measured. For more details see Table 3. 



companies the presence of abundant warm molecular gas 
(Boogert et al. 2000a; van dcr Tak et al. 2000). The re- 
sult obtained for Elias 29 can only be understood once the 
location of the ices and the physical conditions of the var- 
ious gas components are known. Therefore, in this paper 
wc try to identify any foreground material, the presence 
of a circumstellar envelope as well as the presence and 
orientation of a circumstellar disk, and the column den- 
sity of each component. We will then address the question 
where the ices are located, and what their relation is to the 
young star. This study will, as a consequence, reveal im- 
portant information on the nature and evolutionary stage 
of Elias 29, which has many interesting and unique prop- 
erties (Elias 1978). 

Details of the single dish and interferometer observa- 
tions are presented in §2, and the maps and spectra are 
decomposed and interpreted in §3. The physical conditions 
are determined for the different components along the line 
of sight, among which are two foreground clouds (§4.1), 
a remnant envelope and face-on disk (§4.2), as well as a 
dense ridge from which Elias 29 probably formed (§4.3). 
The gas phase conditions and abundances are linked to 
the ice observations in §4.4. The depletion of gas phase 
species is compared to young class objects and quiescent 
clouds in §4.5. The results are summarized in §5. 

2. OBSERVATIONS 

Single dish and interferometer millimeter wave spec- 
tral line and continuum observations were made toward 
the infrared position of Ehas 29 (a = 16^27™09".5, d = 
-24°37'18"; J2000). Rotational lines of CO, HCO+, CS, 
H2CO, CH3OH, and isotopes were selected in the 70-400 
GHz spectral range based on their sensitivity to column 
density, temperature, and density (Blake et al. 1995). 
Below we discuss the technical details of these observa- 
tions. Not discussed is a map of the CO 6-5 line (692 
GHz), for which wc refer to C. Ccccarclli et al. (in prepa- 
ration). Also, the technical details of the 1.3 mm IRAM- 
30m continuum map at 15" spatial resolution that is used 
in this paper are discussed elsewhere (Motte, Andre, & 
Neri 1998). Finally, the infrared spectral energy distribu- 
tion (SED) that we use consists of 2-45, and 45-200 /im 
spectra that were obtained with the ISO-SWS and LWS 

CSO beam parameters are available at http://www.submm.caltech.i 



instruments in apertures of ^25" and ^80" respectively, 
and were extensively discussed in Boogert et al. (2000b). 
A summary of all the observations used in this paper is 
given in Table 1. 

2.1. Singe Dish 70-400 GHz Spectral Line Observations 

Single dish NRA0-12m, JCMT, and CSO miUimeter 
wave observations were made with a single pointing to- 
ward Elias 29 during a number of observing runs in the 
1995-2001 period (Table 1). In some lines, small maps 
were made of the Elias 29 environment with the CSO. At 
the JCMT and CSO, we used the DAS and AOS spectrom- 
eter back-ends in the highest available spectral resolution 
mode (~0.1 km s~^), in the 200-400 GHz spectral range. 
At the NRA0-12m, we used the autocorrelator back-end, 
with a channel width of 0.049 MHz, resulting in Nyquist 
sampled spectra with resolutions of 0.40-0.20 km s""*^ at 
the frequencies 70-145 GHz. 

For low frequency transitions, that trace extended cloud 
structure, we took a = 16'^23'"0P.5, S = -24°36'58" 
(J2000) as an off position, which is found to have very 
little or no ^^CO emission (Loren 1989). For higher fre- 
quency CO transitions, we took an off position of 2700" in 
azimuth, and for other molecules the azimuth offset was 
900". No contamination by emission in the off position is 
evident in our spectra. 

Our line observations were corrected for atmospheric at- 
tenuation and telescope losses, using the standard chop- 
ping wheel technique (Kutner & Ulich 1981). TheNRAO- 
12m data retrieved from the telescope are on a scale, 
and the CSO and JCMT data in T^. To convert to the 
main beam brightness temperature (T^h), we divide by 
the beam efficiencies: 

Tmb = T^/V*ra (NRAO - 12m) (1) 
n/Vmb (JCMT, CSO). (2) 

This corrects for losses in the beam side lobes. For the 
JCMT and CSO, losses due to forward spillover and scat- 
tering are included in r/mh- The beam efficiencies applied 

to our data are summarized in Table 2. For the CSO 
telescope^", r^mb and the main beam size (Table 3) were 

/cso/receivers/beams.html 
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Table 2 
Applied Beam Efficiencies 



Telescope 


Frequency 

GHz 




NRAO-12m 


70-85 


1.00 




85-90 


0.95 




95-115 


0.93 




140-150 


0.81 


JCMT 


200-270 


0.68 




320-380 


0.60 


cso 


200-270 


0.76'' 




320-380 


0.78*' 



^ ri^ given for NRAO-12m 
for 03/1999; 10 % lower in 01/2001 



dard way in the OVRO/MMA reduction package. The 
nearby source NRAO 530 was used as a gain calibrator, 
and flux calibration was performed on the planets Uranus 
and Neptune. The calibrated tracks were combined, and 
deconvolved images were constructed with the MIRIAD 
software package using uniform weighting. The resulting 
synthesized elliptical beam sizes are 4x8" at 87 GHz and 
3x6" at 110 GHz; the noise level achieved per channel in 
the final data is 0.26 Jy/beam and 0.13 Jy/beam respec- 
tively Strong and extended ^^^CO and HCO+ 1-0 emission 
is detected, with perhaps a weak detection of C^^O 1-0 
(Figs. 4-6). SiO 2-1 is undetected toward Elias 29. 

No obvious strong continuum source was seen in the 
OVRO 2.7 and 3.4 mm maps. A weak 3 a peak (7 mJy) is 
detected at 3.4 mm. This is the strongest peak in the map, 
and it is centered on the infrared position of Elias 29. At 
2.7 mm a peak of only 2 cr (5 mJy) significance appears 
at the same position. For our analysis, we take a weighted 
mean of these values, which we will refer to as the con- 
tinuum flux at 3.0 mm: _F(3 mm)=6.1±1.7 mJy. This is 
only a 3.5 a result, and needs confirmation by additional 
observations. 

3. RESULTS 

3.1. Large Scale Emission: Single Dish Spectra 

The single dish maps show that Elias 29 is not a particu- 
larly prominent center of molecular line emission. Within 
a radius of ~ 1.5', the HCO+ 3-2 and C^^O 2-1 emission is 
highly structured, and quite differently distributed (Figs. 
1 and 2). 

The HCO+ 3-2 emission, tracing high densities (> 
10^ cm~^), is concentrated in a remarkable ridge-like 
structure, oriented in the southeast-northwest direction 
(Fig. 1), at a velocity of Vlsr~5.0 km s'^ (Fig. 2). It is 
likely no coincidence that the protostars Elias 29, WL 20, 
LFAM 26, GY 210 as well as the 1.3 mm continuum proto- 
stellar condensations E-MM3 and E-MM5 are all located 
along this dense ridge (Fig. 1). The ridge is also particu- 
larly prominent in the 800 firsi continuum (Johnstone et al. 
2000b). Star formation along dense filamentary structures 
is common in the p Oph cloud, and has been explained 
by the presence of magnetic field tubes or, more likely, by 
externally induced shocks (see Motte et al. 1998 for a 

User's Manual for the NRAO-12m Millimeter- Wave Telescope (Mangum 1999) is available at http://www.tuc.nrao.edu/12meter/obsinfo.html 
JCMT beam parameters are available at http://www.jax;h.hawaii.edu/JACpublic/JCMT/ 



determined by observing Mars during observing runs in 
March 1999 and January 2001. For the NRA0-12m", we 
used the theoretical values of rj^ from the Ruze equation, 
increased by a small factor (~ 1.08), needed to reproduce 
the experimental values (Mangum 1999). As a check, we 
determined rj^ from our observation of Jupiter, yielding 
ry^=0.80 at 145 GHz, in excellent agreement with the as- 
sumed value. For the JCMT, the efficiency factors are 
taken from regularly measured values reported on the in- 
ternets^. 

The calibration accuracy between different observing 
runs is not expected to be better than 25% (e.g. Mangum 
1993). As a check on the calibration accuracy, we ob- 
served the very bright nearby source IRAS 16293-2422 
during our JCMT and CSO runs. We find that line inten- 
sities, observed with the same telescope during the night 
and in separate observing runs, indeed scatter with 20- 
30% variations around the mean. Similar differences are 
found when comparing JCMT observations to the stan- 
dard spectra available at the JCMT internet page. In this 
paper, we will assume an uncertainty of 25% in the line 
intensity, unless noted otherwise. 

The data were reduced with the GILDAS/CLASS re- 
duction package, applying low order baselines and us- 
ing weighted means in averaging individual spectra. The 
reduced single dish spectra and maps are presented in 
Figs. 1- 3 and Table 3, and analyzed in §3.1. 

2.2. Interferometer Spectral Line and Continuum 
Observations 

Elias 29 was observed with the OVRO millimeter array 
in the 1999/2000 season (Table 1). The low elevation of 
the source allowed for only short tracks (^7 hours). The 
digital correlator system was used in the "1 8 2" mode, 
providing simultaneous observations of the C^^O 1-0 and 
""^■^CO 1-0 lines in one local oscillator setting as well as 
the HCO"*" 1-0 and SiO 2-1 lines in another. The spectral 
resolution is 0.35 km s~^, with a velocity coverage of 20 
km . The telescope configurations L, H, and E were 
used in all lines, as well as the continuum. The weather 
conditions were best during the L and H observations, re- 
sulting in Tsys - 500 K at 110 GHz, and 800 K at 87 GHz. 
System temperatures were worse in the E track: 700 and 
1300 K respectively. The data were reduced in a stan- 
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ARA (arcsec) 

Fig. 1. — Overlay of IRAM-30m 1.3 mm continuum map (gray scale and two white contours; Motte et al. 1998) and integrated CSO 
HCO+ 3-2 emission (black contours) of the Elias 29 region, with source identifications. Contour levels for the continuum map are at the 7, 
and 14 c confidence level, with <t = 10 mJy/15" beam, and / Tmb dv = 1.1, 1.6, 4.6 K.km s"! for HCO+ 3-2 (cr = 0.2 K.km s"!). The 
1.3 mm map was aligned to the radio continuum positions of Elias 29 and WL 20. The HCO+ map in turn was aligned to the 1.3 mm map, 
for which a small shift of 5" to the southwest was needed. 



short discussion). 

The C^^O 2-1 emission, a column rather than volume 

density tracer, shows that at least three clouds are present 
along the line of sight of Elias 29 (Fig. 2). The channel 
maps show a cloud at Vlsr^ 2.7 km s~^ that peaks to the 
northeast of Elias 29, and a cloud at Vlsr~ 3.8 km s~^ 
spread rather evenly over the map. The brightest cloud at 
'^5 km peaks prominently near the south-southwest of 
the map, and is probably that in which the dense HCO"*" 



ridge resides, given the similar velocities. All these clouds 
are likely present in the foreground, since absorption in 
the ^^CO emission lines is seen at these velocities (Fig. 3). 

In various other single dish lines these clouds show up 
as discrete emission components (Fig. 3). The 2.7 and 
3.8 km s~^ components are primarily seen in the low- 
lying transitions of the CS, HCO+, H2CO, and CH3OH 
molecules, indicating low densities and temperatures of 
these extended clouds. The emission at 5.0 km s^^, iden- 
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Fig. 2. — Channel maps of C^^^O 2-1 (a) and HCO+ 3-2 (b) emission in the Elias 29 region, observed with the CSO telescope. The emission 
is integrated over a 1 kin s^^ interval, with the central velocity (in km s^'^) indicated in the upper left corners. The contour levels are equally 
spaced / Tmb dv = 0.8, 2.1, ...9.9 K.km s-^ for G^^O (a = 0.20 K.km s-^), and 0.3, 0.9, 2.7 K.km s"! for HCO+ (a = 0.09 K.km s-i). 
Note that the map scale is 2x2 arcmin for C^*0, and 3x3 arcmin for HCO+. Foreground clouds are seen in C^^O 2-1 at 2.5 and 3.7 km s~^. 
The Elias 29 core and surrounding ridge are seen at higher velocities, especially in HCO+. 
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Fig. 3. — Emission lines of CO, CS, HCO+, H2CO, and CH3OH toward Elias 29. Vertical dashed lines indicate diflferent dynamic compo- 
nents at A4jSR=2.7, 3.8 and 5.0 km s~^. The spectra have been arbitrarily scaled and shifted along the intensity axis. See Table 3 for line 
intensities. 
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Table 3 

Rotational Transitions Observed toward Elias 29 (Center Position) with Single Dish Telescopes 



Molecule 


Transition 


Frequency 




f T' J, .a 

J J^MB av 


r WrlM 


I'LSR 


Beam 


Telescope 


Date 






MHz 




K km 


km s 


km s 


arcsec 






CO 


2-1 


230538.0 


15.9 


19.4 


6 


2/6.5 


21 


JCMT 


03/1997 




3-2 


345796.0 


26 


98 


7 


2/6.5 


14 


JCMT 


02/1996 




6-5 




9 


46 


6 


2/6.5 


21 


CSO 


07/2000 




691473.0 


17 


124 


12 


1.8/6.5 


7 


JUMl 


04/1995 


13 CO 


6-5^ 


661067.4 


10 


39 


3.6 


4.78 


7 


JCMT 


04/1995 


Ci^O 


1-0 triplet 


112358.7 


0.40 


0.99 


2.3 




56 


NRAO 


05/1995 






112359.0 


0.94 


1.00 


1.0 


3.58 




NRAO 


05/1995 






112360.0 


0.68 


1.05 


1.5 






NRAO 


05/1995 




2-1 multiplet 


224714.3 


1.76 


4.95 


2.63 


4.08 


22 


JCMT 


03/1995 


c^^o 


1-0 


109782.2 


6.34 


11.0 


2.0 


3.6 


57 


NRAO 


05/1995 




2-1 


219560.4 


5.81 


16.2 


2.62 


4.19 


22 


JCMT 


03/1995 








9.6 


21.3 


2.41 


4.05 


35 


CSO 


01/2001 




3-2 


329330.6 


4.3 


10.4 


2.28 


4.24 


15 


JCMT 


02/1996 


cs 


2-1 


97981.0 


1.28 


3.55 


2.61 


3.96 


64 


NRAO 


05/1995 




5-4 


244935.6 


0.56 


0.94 


1.56 


4.92 


20 


JCMT 


03/1995 




7-6 


342883.0 


<0.09 


< 0.13 






14 


JCMT 


02/1996 


C34S 


2-1 


96412.9 


0.14 


0.23 


1.62 


3.29 


65 


NRAO 


05/1995 


H2CO para 


loi — Ooo 


72838.0 


0.62 


1.44 


2.19 


3.80/5.39 


86 


NRAO 


05/1995 




2o2 — loi 


145603.0 


0.56 


0.89 


2.0 


3.67/4.10 


43 


NRAO 


05/1995 




303 — 2o2 


218222.2 


0.39 


0.66 


1.89 


4.60 


22 


JCMT 


03/1995 




322 — 221 


218475.6 


<0.06 








22 


JCMT 


03/1995 


H2CO ortho 


2l2 — 111 


140839.5 


0.80 


1.85 


2.17 


4.20 


45 


NRAO 


05/1995 




2ll — lio 


150498.4 


0.80 


1.85 


2.17 


4.20 


42 


NRAO 


05/1995 




3i2 — 2ii 


225697.8 


0.41 


0.76 


1.76 


4.93 


21 


JCMT 


03/1995 








0.26 


0.48 


1.70 


4.51 


32 


CSO 


07/2000 


CH3OH 


2-1 triplet 


96739.4 


0.15 


0.19 


1.24 




65 


NRAO 


05/1995 






96741.4 


0.18 


0.33 


1.68 


3.59 




NRAO 


05/1995 






96744.6 


<0.03 










NRAO 


05/1995 




5-4 multiplet 


241802.0 


<0.05 








20 


JCMT 


03/1997 








<0.02 








31 


CSO 


06/2000 


HCO+ 


1-0 


89188.5 


3.57 


5.69 


3.3/1.0 


4.3/4.7 


70 


NRAO 


05/1995 




3-2 


267557.6 


1.06 


2.7 


2.28 


4.77 


18 


JCMT 


05/1996 








4.34 


4.9 


1.05 


4.63 


28 


CSO 


06/2000 




4-3 


356734.3 


0.30 


0.80 


2.5 


4.5 


21 


CSO 


03/1999 


H13CO+ 


1-0 


86754.3 


0.38 


0.36 


0.90 


4.66 


72 


NRAO 


05/1995 




3-2 


260255.5 


0.04 (0.02) 


0.09 (0.02) 


2.15 


4.89 


29 


CSO 


03/1999 



"Calibration errors are 25%, unless noted otherwise in parentheses. 

''CO 6-5 spectra presented and analyzed in detail in C. Ceccarelli et al., in prep. 
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Table 4 

Integrated Intensities of Decomposed Lines at Center Position 



Molecule Transition J Tmb dv for velocity component'' 

2.7 km s-^ 3.8 km s"^ 5.0 km s'^ 





1-0 




1.78 


5.18 


3.97 




2-1 




1.93 


4.19 


9.50 




3-2 




1.00 


2.75 


6.78 




1-0 




0.66 


1.94 


< 1.6 


cs 


2-1 




0.69 


1.01 


1.72 




5-4 




< 0.07 


< 0.09 


0.82 




7-6 




< 0.13 


< 0.13 


< 0.13 




2-1 




0.09 


0.12 


< 0.05 


H2CO para 


loi - 


-Coo 


< 0.15 


0.55 


0.65 




2o2 - 


-loi 


< 0.07 


0.43 


0.39 




3o3 - 


- 2o2 


< 0.03 


< 0.04 


0.82 




322 - 


-221 


< 0.03 


< 0.04 


< 0.12 


H2CO ortho 


2l2 - 


- Ill 


0.30 


0.50 


1.03 




2ii - 


- lio 


0.21 


0.52 


0.65 




3l2 - 


-2ii 


< 0.06 


< 0.09 


0.78 


CH3OH 


2-1 96739.4 


0.03 


0.14 


< 0.05 




2-1 96741.4 


0.06 


0.18 


0.10 




2-1 96744.6 


< 0.02 


< 0.02 


< 0.05 


HCO+ 


1-0 




0.49 


0.46 


4.26 




3-2 




< 0.04 


< 0.07 


3.17 




4-3 




< 0.05 


< 0.08 


0.76 


H13CO+ 


1-0 




< 0.02 


0.05 


0.33 



* Relative uncertainties are smaller than the error in the absolute calibration 



(25%) 

^ Derived from 112360.0 MHz fine structure line, multiplied by a factor of 3 
to correct for emission in other fine structure lines. 



tified with the EHas 29 core and ridge, shows up in the 
high excitation hnes. 

Given the fact that these emission components have dif- 
ferent spatial distributions, their relative intensity also de- 
pends on beam size. The HCO+ 3-2 line is a factor of 4 
stronger and a factor of 2 narrower in the CSO data com- 
pared to JCMT: the larger CSO beam (28" versus 18"; 
Table 3) picks up more emission from the ridge, which 
is bright in the west (§3.2). The wide beam NRA0-12m 
HCO"'" 1-0 and H^^CO"'" 1-0 spectra have remarkably nar- 
row 5.0 km s~^ components, originating from the extended 
ridge. In contrast, the H2CO 3i2 — 2ii line may peak on 
the Elias 29 core, rather than the ridge. It has the same 
width in the CSO and JCMT beams, but the CSO spec- 
trum is weaker (Table 3). 

We have performed a Gaussian decomposition of the 
emission lines to separate the foreground clouds from 
the dense material around Elias 29. The derived rela- 
tive strength of these blended lines depends sensitively on 
the assumed line width. The most reasonable solution, 
based on the H2CO and CS lines and the C^^O map, is 
to simultaneously fit the peak intensity with three Gaus- 
sians, centered on I^lsr^ 2.7, 3.8, and 5.0 km s^^, with 
FWHM=0.5, 0.8, and 1.5 km s^^. When needed, we al- 
low a 30% variation in the FWHM, and 0.20 km s'^ in 
Vlsr- Some of these small variations may be real, but 
given the complexity of the line profiles and spatial dis- 
tribution of the various components, we will not seek a 
physical interpretation for this. The integrated intensities 
are summarized in Table 4, and the physical conditions 



are derived from these in §4. 

We conclude that when studying individual protostars 
in the p Oph cloud, it must be realized that many different 
physical components are present within single dish beams 
of ~ 15" or larger. The spatial information of interferom- 
eters is essential here. 

3.2. Small Scale Emission: Interferometer Maps 

In the OVRO interferometer maps, bright ^"^CO and 
HCO"*" 1-0 emission is detected in the direct neighbour- 
hood of Elias 29 (Fig. 4). The ^^CO emission is most 
strongly peaked on the infrared position, and with a spa- 
tial FWHM- 6" (900 AU) it is resolved along the EW 
direction, where the OVRO beam is smallest (3"). This is 
the Elias 29 core, which likely consists of a disk/envelope 
system (§4.2). An extension of ~ 12" toward the southwest 
is visible in both lines, but most prominently in HCO"'" 1-0, 
which must be attributed to the ridge discussed above. In- 
deed, on this small scale the HCO"*" emission is also parallel 
to the IRAM-30m 1.3 mm continuum emission (Fig. 4). 

The OVRO ^3 CO and HCO"^ 1-0 spectra peak between 
^LSR=5-7 km s~^, and show no evidence for the strong, 
extended 2.7 and 3.8 km s"'^ foreground emission seen in 
the single dish maps (Fig. 5). The C^^O 1-0 emission 
is spectacularly absent in the interferometer spectrum; it 
is at least a factor of 10 weaker compared to the single 
dish data. After dilution to the NRA0-12m beam (tak- 
ing the emitting area from the OVRO ^^CO image), the 
OVRO C^^O signal is even a factor of 60 weaker than 
the bright detection with the NRA0-12m telescope. This 




30 20 10 -10 -20 -30 30 20 10 -10 -20 -30 

ARA (arcsec) ARA (arcsec) 



Fig. 4. — Integrated OVRO "''^CO 1-0 (a) and HCO+ 1-0 (b) emission in the Elias 29 region (white contours), superposed on the 1.3 mm 
contirmum map (gray scale and black contours). Contour levels are drawn at -3,-2, 2, 3, 6, ...,15 and -3,-2, 2. 3, and 4 times the noise level 
(0.48 and 0.90 Jy/beam for ^'-^CO 1-0 and HCO+ 1-0 resp.). The contours for negative signals are given by white dashed lines. Contour levels 
for the 1.3 mm map axe the same as in Fig. 1. The synthesized beam size is given in the top left corners. 




-5 5 10 15 



Fig. 5. — Average line profiles of the central region with more than 2<t emission in the OVRO maps, which is an area of 380 arcsec'^ for 
^•'CO and 650 arcsec^ for HCO+. For the non-detected C-'^^O and SiO lines we took the central 400 arcsec'^. The brightness temperature 
is correspondingly diluted by factors of 7.4, 5.9, 5.8, and 6.4 from top to bottom respectively to derive the expected signal in the single 
dish NRAO-12m beam. The smooth gray line for C^**0 and HCO+ is the NRAO-12m single dish emission, divided by factors of 13 and 9 
respectively (for clarity). The vertical dashed lines represent the single dish emission line components shown in Fig. 3. 




Fig. 6. — Channel maps of OVRO ^■'CO 1-0 (a) and HCO+ 1-0 (b) emission in the Elias 29 region. The map size is 50". The bullets 
indicate the map centers, and facilitate comparison of the panels. Contour levels are drawn at -2, 2, 3, 6, 9, 12 (^^CO) and -2, 2, 3, 4, and 
5 times the noise level (0.13 and 0.26 Jy/beam rcsp.), with the negative contour indicated with dots. Extended emission toward the west is 
present at velocities near 5.0 km s~^ in both lines and coincides with the large scale dense ridge seen in our HCO+ 3-2 single dish map. 
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shows again that the bulk of the molecular material in the 
Elias 29 line of sight is not associated with Elias 29, but 
instead is present in extended foreground clouds, which 
are resolved out by the OVRO interferometer. 

It is interesting to note that the OVRO line profiles arc 
double peaked (Fig. 5), with a 5.0 km s~^ component 
emitting in the dense ridge (Fig. 6). Emission at the cen- 
tral infrared position is strongest at 6 km s~^. Whether 
this is a true dynamical difference between the Elias 29 
core and the ridge, or artificially created by self-absorption 
at ~ 5.5 km s~^, is difficult to answer at present. In this 
paper we will keep referring to this as the 5.0 km s~^ 
component. A disk/envelope/ridge decomposition is at- 
tempted in §4.2. 

4. DISCUSSION: PHYSICAL CONDITIONS 

It is our aim to derive the physical structure of the sur- 
roundings of Elias 29, and to locate the origin of the ices, 
seen abundantly along this line of sight (Boogert et al. 
2000b). In §4.1 the extended clouds at 2.7 and 3.8 km s"^ 
identified above arc discussed; §§4.2 and 4.3 describe the 
more immediate circumstellar environment of Elias 29 in 
terms of a near face-on disk, a remnant envelope, and the 
dense ridge from which the star may have formed (see 
Fig. 7). The physical conditions in each of the compo- 
nents are constrained using the intensities and intensity 
ratios of the single dish and interferometer line emission, 
the 1.3 mm continuum emission (Motte et al. 1998), and 
the infrared SED (Boogert et al. 2000b). This is linked 
to the ice observations in §4.4. 

4.1. Foreground Clouds 

The C^^O line emission at 2.7 and 3.8 km s~^ is ex- 
tended over several arcminutes (Fig. 2), and thus is asso- 
ciated with the overall p Oph cloud complex. Assuming 
that the clouds are sufficiently homogeneous, we can use 



the ratios of the decomposed line intensities (Table 4) to 
find the density and temperature, and subsequently the 
column density, in each of the clouds. We use the escape- 
probability method described by Jansen (1995) to calcu- 
late the molecular excitation. 

The most useful constraints on the density and tempera- 
ture arc given by the intensity ratios of 1-0/3-2 and 
HCO+ 1-0/3-2. Fig. 8 visualizes the density and temper- 
ature values allowed by the observed ratios, taking into ac- 
count line-opacity effects. The 2.7 km s^^ component has 
a temperature T = 15±5 K, and a density n(H2)=(l±0.5) 
xlO^ cm~^. For the 3.8 km s~^ component, the den- 
sity is less than 10''' cm~'^, but the temperature cannot 
be significantly constrained from the line ratios. Here we 
use the far infrared SED to find that Tkin = 25 ± 15 K 
(§4.2). Taking a C^^O abimdance with respect to H2 of 
3.6 X 10-^ (iV[CO]/iV[Ci*O]=560, Wilson & Rood 1994; 
A^[H2]/N[CO]=5000, Lacy et al. 1994), the H2 column 
densities of the 2.7 and 3.8 km s"^ clouds are 5 x 10^^ 
cm~^ and 1.4 x 10^^ cm~^, respectively (Table 5). This 
also assumes that CO is not strongly depleted in these 
clouds, which is validated in §4.5. These parameters fit 
the optically thin C^^O lines also, indicating that a cor- 
rect line opacity was used 

With these physical parameters at hand, we calculate 
that the optical depth of the clouds in the CO 5-6 transi- 
tion is ~4. This is sufficiently high to explain the absorp- 
tion in the CO 6-5 lines at 2.7 and 3.8 km s^^ (Fig. 3). 
The bright CO 6-5 emission is closely associated with Elias 
29 (C. Ceccarelli et al., in prep.), and therefore the 2.7 and 
3.8 km s""'^ clouds are located in front of this source. 

Independent information on the physical conditions in 
these extended foreground clouds is obtained from the 
SED longward of 100 /xm (§4.2). The columns and tem- 
peratures compare also favorably with the extended cold 
dust found by sensitive 180-1100 fim balloon-based mea- 




FiG. 7. — Schematic overview of the proposed geometry of Elias 29 and its environment (not to scale). The line of sight to Elias 29 may 
cross part of the envelope if the system is slightly more inclined than shown here. Most of the ice-absorption bands originate in the foreground 
clouds. 
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surements (Ristorcelli et al. 1999; Table 6). 

The observed lines of CS, H2CO, and CH3OH do not 
further constrain the physical conditions, but they can be 
used to derive abundances. These are listed in Table 5, 
and further discussed in §4.5. 
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Fig. 8. — Temperature and density range constrained by intensity 
ratios of Ci*0 1-0/3-2 (thin lines) and HCO+ 1-0/3-2 (thick lines) 
for the foreground clouds at 1^SR=2.7 (a) and 3.8 km (b), and 
for the ridge material at 5.0 km s~^ (c). The intensity ratios for 
the ridge material have been corrected for the expected contribution 
from the envelope and disk of Elias 29. The filled star indicates the 
adopted temperature and density for each of the components. 



4.2. The Circumstellar Environment of Elias 29 

Understanding the circumstellar environment of Elias 29 
requires the combination of several pieces of crucial in- 
formation: the single-dish HCO+ 3-2 map and 1.3 mm 
dust-continuum distribution (Fig. 1), the spatially resolved 
emission in ^^CO 1-0 and HCO+ 1-0 observed by OVRO 
(Fig. 4), and the infrared SED obtained by the SWS and 
LWS spectrometers of the ISO satellite (Boogert et al. 
2000b; Fig. 9). 

The HCO+ 3-2 and 1.3 mm continuum maps show that 
Elias 29 is located in a narrow, dense ridge, < 30" wide 
and several arcminutes long. The strongest HCO+ 3 -2 
peak in the ridge is situated near Elias 29. The OVRO 
images and channel maps (Figs. 4 and 6) show that part 
of the ^^CO and HCO"'" emission is centered on Elias 29, 
while a significant fraction follows the crest of the ridge, 
10"-20" offset to the west /southwest. This shows that 
multiple components are present in the immediate vicin- 
ity of Elias 29 (i.e., within a typical single-dish beam), 
in addition to the multiple velocity components along the 
line of sight (§4.1). The velocity coincidence between the 
ridge and Elias 29 further supports that the star formed 
from the ridge. 

4.2.1. SED and ConfAnuum, Modeling 

The most important clues about the nature of Elias 29 
are offered by the infrared SED. Because of the large beam 



with which these measurements were taken (§2), the SED 
reflects not only Elias 29 but also many of the other com- 
ponents identified above. Our ability to derive accurate 
parameters for each component therefore proved essential 
in helping to understand the nature of Elias 29. The SED 
(Fig. 9) is remarkably flat between 10 and 200 /xm. At 10 
yum there is a prominent silicate absorption band, and at 
5 jim the SED shows a broad peak on top of which a num- 
ber of sharp ice absorption features are present. Below 2 
/im the emission drops off sharply. Flat protostellar SEDs 
are generally explained by a circumstellar disk. Compared 
to other circumstellar material distributions, a disk, espe- 
cially a flared one, has the surface-area vs. temperature 
distribution required to create a flat SED. 

For our analysis, we adopt the flared disk model of Chi- 
ang & Goldreich (1997) with the parameters listed in Ta- 
ble 7. The flatness of the SED is explained in this model 
by a superheated surface layer. The continuum emission of 
the disk is modeled using the radiative transfer part of the 
Monte-Carlo code of Hogerheijde & van der Tak (2000); 
this part of the code simply calculates the expected con- 
tinuum emission from a given distribution of density and 
temperature by building up a grid of points on the sky, 
for each of which the radiative transfer is solved along 
straight lines. The resulting grid is then convolved with 
the appropriate beam sizes. Because of the large number 
of parameters involved, their values should be considered 
as representative only. We believe that the general charac- 
teristics of the model are firmly established, however. The 
flatness of the SED over a large wavelength range (4-100 
^m; Fig. 9) limits the inclination of the system to < 60° 
(90° being edge-on; Chiang & Goldreich 1999). Such a low 
inclination cannot explain the deep ice and silicate absorp- 
tion features that are visible in the observed SED; these 
must originate in the envelope and foreground clouds. The 
low disk inclination is compatible with the very high veloc- 
ity 80 km s^^), variable outflowing hot CO gas seen in 
this object (A.C.A. Boogert, G.A. Blake, & M.R. Hoger- 
heijde, in prep.). A fully face on orientation is however not 
likely, given the presence of low surface brightness scat- 
tered K band hght out to distances of 15" (2400 AU) from 
the central object (Zinnecker, Perrier, & Chelli 1988). This 
light presumably traces the outflow lobes, adjacent to the 
remnant envelope and the outer edges of the disk. The 
emission is extended along the southwest/northeast direc- 
tion, perpendicular to the high velocity outflowing gas (C. 
Ceccarelli et al., in prep.), with an axis ratio compatible 
with an inclination of roughly ~ 30°. Flattening on a sim- 
ilar scale is observed in the OVRO images and in the 1.3 
mm continuum IRAM-30m single dish map. 

The flux between 20 and 40 /im indicates a tempera- 
ture scaled upward by a factor of 2.25 with respect to the 
values used by Chiang & Goldreich (1997). The central 
star is therefore approximately 50 times brighter than their 
standard model, presumably because it has a higher mass 
(using the scaling relation T oc L^/'"'; Chandler & Richer 
2000) . The disk midplane dominates the emission beyond 
25 /txm, and the shape of the SED at these wavelengths 
limits the disk radius to ^ 500 AU. The warmer surface 
layer dominates the shorter wavelengths and generate a sil- 
icate emission feature at 10 /zm (Fig. 9). The disk model 
can explain only ~20% of the emission maximum around 
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Table 5 

Physical Conditions in Elias 29 Foreground Clouds and Ridge, and Comparison with Other Clouds 



Quantity Unit Cloud component TMC I'' NGC 1333' 







2.7 km 


3.8 km 


5.0 km 




IRAS 4A 


IRAS 4B 


bo 


[km s"^] 


0.3 


0.5 


1.0 








Tkin 


[K] 

[10'' cm-3] 


15 (5) 


25 (15) 


15 (5) 








n(H2) 


1.0 (0.5) 


< 10 


40 (20) 








iV(H2) 


[10^2 cm-2] 


0.5 (0.2) 


1.4 (0.2) 


1.0 (0.2) 


1.0 


14 


6 


X(C'®0)" 


[10-9] 


360 


360 


360 


304 


40 


70 


X(Ci^O)" 


[10-9] 


110 


110 


110 


95 


7.6 


19 


X{CS) 


[10-9] 


6.0 


15.0 


6 


6 


1.2 


0.2 


X(HCO+) 


[10-9] 


0.6 


0.7 


b 


9 


0.4 


0.1 


X(ortho H2CO) 


[10-9] 


3.0 


8.0 


b 


50 


0.4 


0.9 


X(para H2CO) 


[10-9] 


<1.0 


2.7 


b 








X(CH30H) 


[10-9] 


0.4 


0.5 


b 


3 






Av 


[mag] 


2.9 


8.2 


< 5.9" 









General: all abundances X are w.r.t. H2 and have an error of 50% 

Prom H2 using standard relations (see text), except TMC-1 
^ Line intensities in ridge have to be corrected for expected contribution from envelope and disk (see text). 
No independent determination of abundances is therefore possible. 

Upper limit; only a fraction of the ridge material may actually be in front of Elias 29. 

on 'CP' peak; Pratap et al. 1997; Ohishi & Kaifu 1998 
^ Blake et al. 1995, assuming iV(H2) from dust 



Table 6 

Column Density of Cold Gas Derived by Different Methods 



Method 


Component 


Ar(H2)=' 

10^2 cm-2 




silicate absorption'^ 




2.5-6 


14-34 


4.7 /nm i^CO absorption'^ 


feD=l-0 km s^-' 


3 (0.5) 


17 (3) 


C^^O 2-1 emission 


6d=0-7 km s^-' 


11 (4) 


63 (23) 


VtsR=2.7 


0.5 (0.2) 


2.9 (1.1) 




VtsR=3.8 


1.4 (0.2) 


8.0 (1.1) 




VtsR=5.0 


1.0 (0.2) 


5.7 (1.1) 


ISO SED and IRAM-30m 1.3 


mm envelope 


0.3-1 


1.7-6 


0.18-1.1 mm continuum'^ 


10-15 K 


0.5-5 


2.9-29 


X rays^ 




3.8 


22 


assuming iV(H2)/Ar(CO) = 


5 X 10"*; conversion 


factor j4v = 


Ar(H2) X 8.6/15 X 10^^ 



Bohlin, Savage, & Drake 1978; '^ see Boogert et al. 2000b; '^ A.C.A. Boogert et al., in 
preparation; " Balloon experiment: Ristorcelli et al. 1999; ^ Imanishi, Koyama, & Tsuboi 
2001 
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A (jU.m) 

Fig. 9. — Infrared SED of Elias 29 as observed by ISO LWS and SWS (thin line with ice and silicate absorption features; Boogert et al. 
2000b) combined with the 1.3 mm flux in the 15" IRAM-30m beam (open square; Motte et al. 1998). The thick gray line shows the modeled 
SED as described in §4.2. It consists of four components: (1) a 0.012 Mq flared circumstellar disk following Chiang & Goldreich (1997) with 
a two-layer vertical temperature structure (cool interior [long-dashed line] and superheated surface [dotted line]); (2) a 0.12 M© collapsing 
envelope (dash-dot line); (3) an ad hoc 800 K central blackbody with a 1.4 AU radius (dash-triple dot line), required to fit the fiux maximum 
at A ~ 5 nm and probably corresponding to the warmest regions of the disk close to the star and/or the stellar photosphere extincted by the 
disk; and (4) a 19 K foreground layer (thin solid line pealcing at ~120 /im) with a column density of iV(H2) = 1.9 x 10^^ cm~-^ as derived for 
the material at 1^SR=2.7 and 3.8 km s~^ in §4.1, which is responsible for the continuum absorption below 3 /um in the final four component 
model, and for the prominent silicate and ice absorption bands. 



5 fxm.. Indeed, speckle observations of Elias 29 reveal the 
presence of a -^400 AU radius thermally emitting region 
(T a few 100 K; i.e. not scattered light) responsible for 
~20% of the M band flux (Zinncckcr ct al. 1988), which, 
in our model, is explained by the warm disk surface layer. 
The observed remaining 80% of the M band emission orig- 
inates from hot dust within a few AU from the central ob- 
ject. Perhaps this is a puffed up hot rim at the inner disk 
edge, just outside the region where dust has evaporated 
(Tsubi ~ 1500 K; Dullemond, Dominik, & Natta 2001; 
Zinnecker et al. 1988). The innermost, hottest, part of 
the rim may be related to the 0.25 AU radius structure 
found in lunar occultation K band observations (Simon ct 
al. 1987). Wc do not attempt to fit the inner regions 
and stellar photosphere with that level of detail, but in- 
stead use an 800 K blackbody with an effective radius of 
1.4 AU to explain the remaining 80% of M band emission, 
and roughly the shape of the SED at shorter wavelengths 
(Fig. 9). 

The disk alone is also insufficient to explain the emis- 
sion beyond ~ 55 /xm. The single-dish line emission clearly 

reveals the presence of appreciable columns of material at 
2.7 and 3.8 km s^^ (§4.1). When filling the ISO beam with 
these columns, the emission between 55 and 180 /xm can 
be fit at a dust temperature of 19 K. This material pro- 



vides sufficient opacity in the 10 ^m silicate band to turn 
the disk's emission feature into the observed absorption 
feature. The opacity shortward of 2 /xm is also sufficient 
to explain the observed steep drop in emission at short 
wavelengths. This strongly suggests that a large fraction 
of the 2.7 and 3.8 km s~^ material is located in front of 
Elias 29, as was found from the absorption in the ^^CO 
emission lines as well (§4.1). 
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Fig. 10. — Radial intensity profile along the ridge centered on 
Elias 29 obtained from the 1.3 mm map of Motte et al. (1998) 
(bullets with error bars), compared to the predicted profile based 
on a model for Elias 29 with a disk and an envelope (solid line) 
and with a disk only (dashed line). Both models include a 19 K, 
Ar(H2) = 6 X 10^1 cm-2 background to fit the levelling off of the 
emission beyond a radius of 25 — 50" . Note that the true total back- 
ground column, filtered out of the IRAM-30m 1.3 mm data, is a 
factor of 3 higher (§4.2.1). 

While the model of disk and foreground clouds fits the 
SED well, it produces a 1.3 mm continuum emission that 
is a factor of 3 larger than observed with the IRAM-30m 
telescope, at large distances from the object 30 mJy; 
Fig. 10). This may result directly from the adopted 'dual- 
beam' mapping procedure of the IRAM-30m observations, 
which filters out emission on the chopping scale (Motte et 
al. 1998; Johnstone et al. 2000a). Additionally, the 
model of disk and foreground predicts a source size at 
1.3 mm that is too small (Fig. 10). The observed pro- 
file indicates that the source extends to a radius of at 
least 3800 AU (25"), much larger than is realistic for a 
circumstellar disk. Instead, it suggests that Elias 29 and 
its disk are embedded in a residual cloud condensation 



from which the system formed. For lack of stronger con- 
straints, we model this envelope (as we will refer to this 
condensation from now on) with the inside-out collapse 
model of Shu (1977). Earlier work (e.g., Hogerheijde & 
Sandell 2000; Ceccarelh et al. 2000) has shown that this 
model provides an adequate description of the density and 
temperature of envelopes around YSOs. The two parame- 
ters of this model are the sound speed and the time since 
the collapse started, for which we take, rather arbitrarily, 
0.13 km s^^ and 1.2 x 10"'' yr respectively. This gives an 
envelope mass of 0.12 Mq and a radius of the 'collapse 
expansion wave' of 6000 AU. Additional parameters are 
listed in Table 7. This model, combined with the disk, 
produces a radial emission profile that fits the measured 
1.3 mm contimmm distribution much better than the disk 
alone (Fig. 10). The SED below ~ 100 /xm is remarkably 
little affected by the envelope, even though with 0.12 
it is much more massive than the disk (0.012 M©; Table 7). 
The material in the disk, of course, has a somewhat higher 
temperature range than the envelope (200-33 K vs. 44-17 
K, respectively), yielding significantly more infrared fltix. 
The envelope produces only a small amount of additional 
absorption in the 10 /im silicate band. 

The parameters for the disk and envelope derived above 
(Table 7) maximize the contribution to the SED from the 
disk and minimize the mass of the envelope. There are 
several pieces of evidence that favor, but not proof, this 
maximum disk model above models with smaller disks. 
First, as for sources with directly detected disks (Hoger- 
heijde 1998), the disk contributes about half of the 1.3 
mm continuum fiux (^^ 120 mJy/15" beam; foregroimd 
subtracted). Indeed, assuming that the weak 3 mm con- 
tinuum fiux within the small OVRO beam (§2.2) orginates 
from the disk, then the expected 1.3 mm disk flux would 
be half the single dish flux (i^[1.3 mm]=49±14 mJy, tak- 
ing spectral index 2.5; Hogerheijde 1998). Second, the 
abovementioned 400 AU radius structure, emitting ther- 
mally at 5 /xm, may well be the surface of a large, warm 
disk. And third, a large face-on flared disk is the simplest 



Table 7 

Parameters of Disk and Envelope Model 



Item 



Value 



Disk: Chiang & Goldreich 1997 



Temperature 
Dust opacity 

Inner radius 
Outer radius 
Mass 

% mass in superheated surface layer 



2.25 X standard model of Chiang & Goldreich 

Ossenkopf & Henning 1994, MRN model with 10^ years of coagulation at 

10^ cm-3 

0.01 AU 

500 AU 

0.012 M0 

2.5% 



Envelope: Shu 1977 



Sound speed 

Age 

Radius of collapse expansion wave 

Temperature 

Dust opacity 

Inner radius 
Outer radius 
Mass 



0.13 km s^i 
2.2 X 10"> yr 
6000 AU 

35(r/1000AU)-0-4 K 

Ossenkopf & Henning 1994, MRN model with 10^ years of coagulation at 

10^ cm~^ 
500 AU (3") 
6000 AU (39") 
0.12 M0 
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way to explain the observed flat infrared SED. 

Nevertheless, the evidence for a large 500 AU radius 
disk is not conclusive and therefore we also try to max- 
imize the envelope's contribution in our model. We find 
that its mass can be as large as 0.33 Mq, while the cor- 
responding disk mass needs to be reduced to 0.002 M©, 
in part by reducing the disk's radius to 30 AU to prefer- 
entially decrease the amount of cool material. The disk's 
temperature scaling does not change appreciably with re- 
spect to that used before. It must be emphasized that in 
this minimum disk scenario, the flatness of the infrared 
SED is the result of a "disk/envelope conspiracy", rather 
than a direct consequence of the flared disk's properties. 
Although the disk and envelope masses arc different from 
those listed in Table 7, they are comparable to within the 
same order of magnitude and this level of uncertainty does 
not affect our conclusions. In the remainder of this section 
wc will adopt the slightly preferred values from Table 7, 
i.e. the maximum disk case. 

4.2.2. Emission Line Modeling 

We use the accelerated Monte- Carlo method of Hoger- 
heijde & van der Tak (2000) to calculate the molecular 
excitation and the resulting line emission in the various 
beams, including a realistic sampling of spatial scales re- 
covered by the interferometer observations. Because the 
disk and envelope exhibit a range of densities and temper- 
atures, we cannot use the much simpler escape probability 
approach here. We calculate the line intensities for three 
different cases (Table 8): for the disk alone, for the enve- 
lope alone, and for the combination of disk-|-envelope tak- 
ing into account optical depth effects. The predicted line 
intensities are meant to be illustrative only. The general 
trends are robust, however. Adopting 'fiducial' molecu- 
lar abundances derived for the two foreground clouds (Ta- 
ble 5), we find that the emission of C^^O observed with 
OVRO is overestimated by at least a factor of 5 by the 
disk-|-envelope model. This can be accounted for by adopt- 
ing a depletion factor of 5-10 for CO in the disk. A large 
CO depletion factor has been seen in many other disks 
as well (Dutrey, Guilloteau, & Guelin 1997; Shuping et 
al. 2001). From the single-dish lines it can also be con- 
cluded, in particular from CS 7-6, that a significant de- 
pletion in the disk is needed (Table 5). Furthermore, we 
find that the higher excitation transitions are dominated 
by the disk-l-envelope, while the lower excitation lines can 
only be partially explained by this model. A fair fraction 
of these lines must therefore originate in the ridge. 

4.3. The Dense Ridge 

As mentioned above, Elias 29 appears to be embedded 
in a ridge traced by the 1.3 mm continuum and HCO"*" 
3-2 line emission, but shghtly offset from the crest of the 
ridge. How much of the line emission measured in the 
single-dish and interferometer beams at 5.0 km s^^ can 
be attributed to the Elias 29 disk-henvelope system, and 
how much originates in this dense ridge? This question 
is relevant, because it addresses the issue of how much 
additional column density resides in the ridge that can 
harbor the ice-absorption bands. When we subtract the 
predicted emission of Elias 29's disk and envelope from the 
single-dish line fluxes at 5.0 km s~^ of Table 4, we can use 



the ratios of the remaining intensities to roughly constrain 
the conditions in the ridge. Fig. 8c shows the results for 
the C^^O 1-0/3-2 and HCO+ 1-0/3-2 lines. The density 
is constrained at (4 ± 2) x 10^ cm~^, significantly higher 
than that of the foreground clouds, while the temperature 
is found to be similarly low at 15 ± 5 K. These param- 
eters are consistent with the appearance of the ridge in 
the 1.3 mm and HCO+ 3-2 maps (cold and dense). The 
corresponding column density along the ridge is 1 x 10^^ 
crn^^, comparable to that of the largest column density 
foreground cloud at 3.8 km s~^. Elias 29 is offset from 
the maximum column density in the ridge, however, and 
the pencil beam traced in the infrared absorption features 
may contain a significantly smaller column density. The 
fact that the foreground clouds and, to a lesser extent, 
the envelope can already fit the 10 /xm absorption would 
suggest that the ridge does not contribute significantly. 
This conclusion does not change if we assume the higher 
envelope mass of 0.3 Mq. 

AA. Location and Thermal History of Lees toward 

Elias 29 

We have identified four regions of significant column 
density in the line of sight of Elias 29: the 2.7 and 3.8 
km s~^ foreground clouds, the dense ridge, and the en- 
velope. All of these are cold environments, in which icy 
mantles can exist, and together they must account for the 
rich spectrum of ices seen with the ISO satellite (Boogert 
et al. 2000b). Their relative contribution to the visual 
extinction and ice bands depends on which fraction of the 
column is actually in front of Elias 29. The 2.7 and 3.8 
km s^^ clouds are likely in front (§4.1) with columns cor- 
responding to visual extinctions ^v=2.9 and 8.2. The 
visual extinction of the dense ridge, Ay <5.9, is an upper 
limit beause the ridge is somewhat displaced from Elias 
29. For the envelope. Ay is 2.5 magn. for the minimum, 
and 6 magn. for the maximum envelope case. These are 
upper limits because the system is seen at low inclination 
and much of the envelope material is likely not in front 
of Elias 29. The extinction from either ridge or envelope 
is however not negligible, given the absorption seen at 5.0 
km s~^ in the CO 6-5 emission line (Fig. 3). Finally, ex- 
tinction by the circumstellar disk is not significant due to 
its low inclination. 

There is observational evidence that ices in molecular 
clouds form only above a certain threshold extinction value 
(^th; Whittet et al. 1988; Tanaka et al. 1990), due 
to photodesorption at the cloud edges (Tielens & Hagen 
1982). Above Ah, the ice abundances grow linearly with 
^v- The visual extinction can thus, in principle, be used 
as an indicator for II2O ice growth and destructions mech- 
anisms. Unfortunately, the usefulness of this relation is 
limited by the uncertainty in Ath- Measurements of a va- 
riety of sight-lines in p Oph indicate Ath=12 (Teixeira & 
Emerson 1999b), but ^th < 7 toward OB stars tracing 
ices in foreground clouds only (Shuping et al. 2000). The 
Ay values used to calibrate the aforementioned relation 
are highly uncertain as well (Teixeira & Emerson 1999b). 
While our emission line method has its own uncertainties 
(beam dilution, disk/envelope ratio), the generally used 
H-K color method may be uncertain by a factor of 2, for 
example due to an unknown color of the intrinsic light 
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Table 8 
Predicted Line Intensities 



Line 


Observed* 


Envelope+Disk'' 


Envelope*^ 


Disk'l 


Adopted abundance 






Single dish observations (K km s 






U yJ i— U 


o.y / 


0.57 


0.52 


n Ac; 
U.UO 


Q R 1 n — 7 
o.D X lU 


{^1S{~\ o 1 
C L> z— 1 


9.50 


3.22 


2.64 


0.75 




U O 3-2 


6.78 


4.57 


2.97 


2.15 




HCO+ 1-0 


4.26 


1.01 


0.99 


0.10 


3 X lO"** 


HCO^ 3-2 


3.17 


2.32 


1.98 


3.43 




HCO+ 4-3 


0.76 


1.51 


0.69 


2.87 




Hi3C0+ 1-0 


0.33 


0.04 


0.03 


0.01 


4 X 10"" 


CS 2-1 


1.72 


0.78 


0.76 


0.13 


1 X 10~® 


CS 5-4 


0.82 


1.49 


0.40 


2.14 




CS 7-6 


< 0.13 


4.44 


0.11 


4.39 




C^^S 2-1 


< 0.05 


0.07 


0.06 


0.02 


4 X 10-1° 


ortho-H2CO 2i2-lii 


1.03 


1.16 


1.13 


0.27 


8 X lO-'' 


ortho-H2CO 2ii-li,) 


0.65 


0.86 


0.83 


0.31 




ortho-H2CO 3i2-2ii 


0.78 


1.0 


0.77 


1.66 




para-H2CO loi-Ooo 


0.65 


0.19 


0.18 


0.04 


2 X 10-8 


para-H2CO 2o2-loi 


0.39 


0.06 


0.02 


0.05 




para-H2CO 3o3-2o2 


0.82 


0.90 


0.54 


0.99 




para-H2CO 322-221 


< 0.12 


1.62 


0.09 


1.55 





Interferometer observations (Jy beam ^, image peak per channel): 



C^^O 1-0 <0.26 1.22 0.14 1.93 

i^CO 1-0 ~ 1.56 1.77 0.82 4.76 

HCO+ 1-0 1.30 1.11 0.92 3.77 

^ intensity entire 5.0 km s— ^ component, including ridge; taken from Table 4 
^ modeled intensity for disk/cnvclopc system: no depletion assumed 

modeled intensity envelope: infinite depletion disk 

modeled intensity disk; infinite depletion envelope 



source. Indeed, for Elias 29 a much lower Ay < 29 was 
found from J-H color (Th.P. Greene, priv. comm.), com- 
pared to the H-K color (Ay = 47; Wilking & Lada 1983). 
Nevertheless, using the relation in Teixeira & Emerson 
(1999b), and assuming that the threshold extinction for 
H2O formation (Ath=12) applies only once for all clouds 
combined, we find that the total Ay < 23 corresponds to 
A'^(H20) < 11 X 10" cm-2, which is <30% of the observed 
ice column density (Boogert et al. 2000b). In this pic- 
ture, the observed ice column is thus much higher than 
that expected from the total visual extinction along the 
line of sight. Given the abovementioned uncertainties in 
the ice column versus Ay and ^th relation, the most reli- 
able conclusion that can be drawn from this is that H2O ice 
sublimation does not play an important role in the Elias 
29 line of sight. 

Prom a different perspective, the low temperature of the 
clouds along the line of sight (Table 5) indicates that in- 
deed no sublimation or other forms of thermal process- 
ing (e.g. crystalhzation) of polar, H20-rich ices can take 
place. A similar conclusion was previously drawn from 
the (low) far infrared color index, which is found to be a 
good indicator of thermal processing of ices (Boogert ct 
al. 2000a; van der Tak et al. 2000). Even the tem- 
perature of the warmest absorbing component, the inner 
part of the envelope (^45 K), is less than the H2O crys- 
tallization (^ 60 K) and sublimation temperatures (^ 90 
K). This explains the observed lack of crystallization sig- 
natures in the absorption band profiles of solid H2O and 
CO2 (Boogert et al. 2000b). 

Thermal processing must have played a role for apolar, 



CO-rich ices toward Elias 29, however. The solid CO/H2O 
ratio is 5%, which is a factor of five less than in several 
other lines of sight within the p Oph cloud (Shuping et 
al. 2000). We find that Elias 29 has heated its enve- 
lope sufficiently (17 44 K; §4.2) to evaporate CO rich ices 
(Tsubi ~15-20 k; Tielens et al. 1991), but preserve H2O- 
rich ices. The envelope, however, harbors only a small 
fraction of the ices, and it is unlikely that this 38 Lq ob- 
ject can heat ices located beyond the envelope, such as 
the dense ridge and the two foreground clouds. Neverthe- 
less, the temperatures of the gas and dust in the quiescent 
foreground clouds are similar to, or slightly higher (Ta- 
ble 5) than the sublimation temperature of apolar ices. 
Several known luminous B stars and the Sco 0B2 asso- 
ciation play an important role in heating the dust in the 
p Oph cloud (Greene & Young 1989). The differences of 
the solid CO/II2O ratio within p Oph may thus be ex- 
plained by the location of the ices with respect to these 
luminous heating sources. Bright external heating sources 
are absent in the Taurus molecular cloud, and indeed the 
CO/H2O ratios are significantly larger in that cloud (Teix- 
eira & Emerson 1999b). 

4.5. Depletion 

The depletion factor of CO (defined as solid/[gas-|-solid]) 
amounts to a few percent, assuming the ices are present 
in several of the clouds toward Elias 29 (using A''[CO- 
ice]=1.7 X 10" cm"2. Boogert et al. 2000b). Direct 
measurements of solid CO yield larger depletion factors 
in the quiescent Taurus molecular cloud (8-40%; Chiar et 
al. 1995). Indirect measurements, that rely on gas phase 
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abundances only, found CO depletion factors as high as 
90% toward deeply embedded class sources such as NGC 
1333: IRAS 4A and 4B (Blake et al. 1995). We will com- 
pare gas phase abundances of other molecules toward Elias 
29 with those of class soiu-ces, as well as the well studied 
dense cloud TMC 1 (Table 5). 

Even when taking the uncertainty of gas phase abun- 
dances into account, we find that CO, CS, H2CO, and 
possibly HCO"*" are an order of magnitude more depleted 
in the class sources compared to the clouds in front of 
Elias 29. This can be ascribed to the much larger den- 
sities and low dust temperatures in the outer regions (> 
700 AU) of these class objects, and thus shorter deple- 
tion times scales (Blake et al. 1995). Depletion is much 
lower in the TMC 1 cloud compared to class objects (Ta- 
ble 5). In fact, the H2CO and CH3OH abundances are an 
order of magnitude larger than toward Elias 29. Whether 
these species are depleted toward Elias 29 cannot be an- 
swered directly, because the upper limit to the ice abun- 
dances are typically a fewxlO^^ (Boogert et al. 2000b; 
taking N{R2) = 2.9 x 10^^ cm'^). This is four and two 
orders of magnitude larger than the gas phase CH3OH and 
II2CO abundances respectively, and demonstrates that a 
direct determination of depletion factors by using infrared 
ice bands suffers from a lack of sensitivity. Alternatively, 
the difference in gas phase CH3OH and H2CO abundances 
between these clouds may have its origin in different chem- 
ical evolutionary states, such was found within the TMC 
1 cloud itself for several carbon bearing species (Pratap 
et al. 1997). Indeed, H2CO abundances of several other 
quiescent and star forming cores (2-10xl0~^; Dickens & 
Irvine 1999) are in much better agreement with Elias 29. 

In some class sources, the CII3OH and H2CO abun- 
dances and temperatures are larger than in our comparison 
source NGC 1333: IRAS 4, possibly because they are more 
evolved (van Dishoeck et al. 1995; Blake et al. 1995). 
This warm gas originates from evaporated ices, that are 
enriched in hydrogenated molecules by grain surface reac- 
tions during the collapse of the cold envelope. Shocks (van 
Dishoeck et al. 1995) and thermal evaporation (Ceccarelli 
et al. 2001) in the inner envelope have been proposed as 
ice desorption mechanisms, which can enhance the H2CO 
abundance at radii < 150 AU by 2 orders of magnitude. 
The single dish observations of Elias 29 do not show ev- 
idence of such warm gas, as is illustrated in Fig. 11 by 
comparing the CII3OII spectrum with the class source 
IRAS 16293-2422. This may well be an effect of the larger 
column and accretion rate of the younger class objects. 
Infrared absorption line studies do reveal, however, that 
significant amounts of warm CO and II2O are present on 
small scales near the protostar, partly related to the out- 
flow, and probably partly present in the accreting disk's 
atmosphere (Boogert et al. 2000b; A.C.A. Boogert et al., 
in preparation). Emission from this gas is severely diluted 
in the large single dish milimeter wave beams, and thus 
difficult to detect. 

Finally, the HCO+ abundance toward Elias 29 is compa- 
rable to that found toward class objects, but an order of 
magnitude lower compared to TMC 1. The importance of 
depletion is more difficult to establish, because the HCO"*" 
abundance is highly dependent on local cosmic ray flux 
and electron density. Even in the direction of Elias 29 



the IICO+ abundance varies dramatically. Although the 
HCO"*" 1-0 and CS 2-1 lines have similar critical densi- 
ties, the ridge (at 5.0 km s~^) is much more prominent in 
HCO+ than in CS, with respect to the 2.7 and 3.8 km s~^ 
foregroimd clouds (Fig. 3). We stress however that the 
prominent appearance of the ridge in the HCO"*" 3-2 map 
(Fig. 1), is both an effect of increased HCO+ abundance 
and increased density: the ridge is weaker in the column 
density tracer C^^O, and stronger in the density tracer CS, 
relative to the foreground clouds. 



IRAS 16293-2482 (Class 0) 




Elias 29 (Class 1) 
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Fig. 11. — Low resolution spectrum showing the enormous differ- 
ence in line intensity between the class object IRAS 16293-2422 
and the class I object Elias 29, in all the lines of SO2, CS, and 
CH3OH. 



5. SUMMARY, CONCLUSIONS AND FUTURE WORK 

We have analyzed infrared and millimeter wave line and 
continuum observations to construct a model of the class 
I protostar Elias 29 and its environment. This model has 
to contain a number of different components (summarized 
in Fig. 7): a disk to account for the ^2 50 /im SED, an 
envelope contributing to the emission at 1.3 mm (in partic- 
ular its size), a dense ridge from which Elias 29 may have 
condensed, and foreground material which provides most 
of the extinction. Elias 29 can then be well described by 
a 500 AU radius face-on flared disk with a mass of 0.012 
Mq, embedded in a 6000 AU radius, 0.12 Mq envelope. 
This large disk provides the simplest explanation for the 
observed flat SED, weakly detected 3 mm continuum emis- 
sion, and 400 AU radius 5 /xm thermal continuum emission. 
The present data does however not fully exclude models 
with smaller disks. The minimum possible disk has a 30 
AU radius and a mass of 0.002 Mq surrounded by an en- 
velope of 0.33 Mq. In this case the combination of disk 
and envelope emission produces a flat SED. 

The entire system is embedded in a long, dense, cold, 
and HCO^-rich ridge. Elias 29 is slightly offset from the 
crest of this ridge. In front of the disk, envelope, and ridge 
system are two foreground clouds at a few km lower 
radial velocities that cover the entire field of view. The 
large column of the foreground clouds, corresponding to 



The low mass protostar Elias 29 



17 



Ay ^11, may be responsible for the 'class F appearance 
of Elias 29, which would otherwise appear as a T Tauri 
or Her big Ae/Be star (i.e., optically visible). These same 
foreground clouds are also the most likely repository of 
most of the ices seen along the line of sight 70%). 
The low temperature of the foreground clouds explains 
the observed absence of crystallized ices, and the presence 
of large abundances of polar, H20-rich ices, i.e. thermal 
processing did not play a major role for the ices toward 
Elias 29. The foreground cloud temperature (25±15 K) 
could, however, be high enough to explain the low abun- 
dance of apolar, volatile CO-rich ices, presumably due to 
the proximity of a number of luminous B type stars. The 
important question as to whether the ices in the disk or 
envelope have experienced thermal processing, as is seen 
in the envelopes of massive objects (Boogert et al. 2000a; 
van der Tak et al. 2000), cannot be addressed given the 
large column of foreground material and the face-on nature 
of the system. 

This work shows the value of spectrally and spatially 
resolved information offered by single-dish and interfero- 
metric molecular gas observations in interpreting infrared 



ISO satellite observations of ices along a pencil beam. It 
shows that, at least for the p Oph cloud, it is crucial to 
disentangle the different physical components along the 
line of sight, since otherwise incorrect conclusions may be 
derived, for example, on the origin and evolution of inter- 
stellar and circumstellar ices. 

Follow up observations with sensitive submillimeter in- 
terferometers (e.g. ALMA) or high frequency (> 400 GHz) 
single dish telescopes with small beams will clearly provide 
essential information on the structure, abundances, and 
depletion factor of species in the disk and outflow of Elias 
29, and its relationship to the envelope and ridge. In turn, 
such studies would provide invaluable information on the 
initial conditions of planet formation. 

We thank Remo Tilanus, Goeran Sandell, and Fred Baas 
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contract AST-9980846. The research of MRH is supported 
by the Miller Institute for Basic Research in Science. 
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